ON THE EQUIVALENCE OF CUBOID EQUATIONS 
AND THEIR FACTOR EQUATIONS. 



RUSLAN SHARIPOV 

Abstract. An Eulor cuboid is a rectangular parallelepiped with integer edges and 
integer face diagonals. An Euler cuboid is called perfect if its space diagonal is also 
integer. Some Euler cuboids are already discovered. As for perfect cuboids, none of 
them is currently known and their non-existence is not yet proved. Euler cuboids and 
perfect cuboids are described by certain systems of Diophantine equations. These 
equations possess an intrinsic S3 symmetry. Recently they were factorized with 
respect to this 53 symmetry and the factor equations were derived. In the present 
paper the factor equations are shown to be equivalent to the original cuboid equations 
regarding the search for perfect cuboids and in selecting Euler cuboids. 



1. Introduction. 

The search for perfect cuboids has the long history. The reader can follow this 
history since 1719 in the references [1 44]. In order to write the cuboidal Diophan- 
tine equations we use the following polynomials: 

po = x\ + x\ + xl - L 2 , p 1 =xj + xl-df, 

2 1 2 i2 2 1 2 j2 \ ' ) 

p 2 = x 3 + x 1 -d 2 , p 3 =x 1 +x 2 -d 3 . 

Here Xi, x 2 , £3 are edges of a cuboid, d\, d 2 , d 3 are its face diagonals, and L is 
its space diagonal. An Euler cuboid is described by a system of three Diophantine 
equations. In terms of the polynomials (1.1) these equations are written as 

Pi = 0, p 2 = 0, P3 =0. (1.2) 

In the case of a perfect cuboid the number of equations is greater by one, i. e. 
instead of the equations (1.2) we write the following system of four equations: 

Po = 0, Pl =0, p 2 = 0, P3 = 0. (1.3) 

The permutation group S3 acts upon the cuboid variables xi, X2, £3, di, d 2 , d 3 , 
and L according to the rules expressed by the formulas 

a(xi) = x all a(di) = d ai , a(L) = L. (1.4) 
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The variables x\, x 2 , X3 and d±, d 2l c?3 are usually arranged into a matrix: 



M = 



Xi 


x 2 


X3 


di 


d 2 


d 3 



(1.5) 



The rules (1.4) means that S3 acts upon the matrix (1.5) by permuting its columns. 
Applying the rules (1.4) to the polynomials (1.1), we derive 

&{Pi)=Pa l , o-(po)=Po- (1.6) 

The polynomials p$, px, p 2 , P3 in (1.1) are treated as elements of the polynomial 
ring x 2 , X3, d\, d 2 , d$, L]. For the sake of brevity we denote 

Q[x u X2,x 3 ,di,d 2 ,d 3 ,L} = q[M,L}, (1.7) 

where M is the matrix given by the formula (1.5). 

Definition 1.1. A polynomial p € Q[M, L) is called multisymmetric if it is invari- 
ant with respect to the action (1.4) of the group S3. 

Multisymmetric polynomials constitute a subring in the ring (1.7). We denote 
this subring through SymQ[M, L}. The formulas (1.6) show that the polynomial po 
belongs to the subring SymQ[A/, L], i. e. it is multisymmetric, while the polynomials 
Pi, P2, P3 are not multisymmetric. Nevertheless, the system of equations (1.2) in 
whole is invariant with respect to the action of the group 53. The same is true for 
the system of equations (1.3). 

The polynomials pi, p 2l P3 generate an ideal in the ring Q[M, L]. It is natural 
to call it the cuboid ideal and denote this ideal through 

Ic = (pi,P2,Pa)- (1-8) 
Similarly, one can define the perfect cuboid ideal 

Ipc = (po,Pi,P2,Ps)- (1-9) 

Each polynomial equation p = with p £ Iq follows from the equations (1.2). 
Therefore such an equation is called a cuboid equation. Similarly, each polynomial 
equation p = with p 6 Ipc follows from the equations (1.3). Such an equation is 
called a perfect cuboid equation. 

The symmetry approach to the equations (1.2) and (1.3) initiated in [45] leads 
to studying the following ideals in the ring of multisymmetric polynomials: 



-^C_sym = J c n SymQ[Af, L], / P c.s ym = Ipc n SymQ[M, L}. (1.10) 

Definition 1.2. A polynomial equation of the form p = with p 6 Ic .sym or with 
p G Ipcsym is called an S3 factor equation for the Euler cuboid equations (1.2) or 
for the perfect cuboid equations (1.3) respectively. 

The ideal lpc_sym from (1.10) was studied in [46] (there it was denoted through 
Isym)- The polynomial po used as a generator in (1.9) is multisymmetric in the 
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sense of the definition f.f . Therefore it belongs to the ideal ipcLsym- In [46] this 
polynomial was denoted through p i : 

pi =p = x\ + x 2 2 +x\ - L 2 . (1.11) 

In addition to (1.11) in [46] the following seven polynomials were considered: 



P2=Pl+P2+P3 = {x\ +xl~ df) + 

+ {xl+xl-dD + ixl+xj-dl), 

P3 = dipi + d 2 p 2 + d 3 p 3 = di (x\ +xj -df) + 
+ d 2 {x\ + x\- di) + d 3 {x\ + x%- d 2 ), 

pi =xxpi + X2P2 + x 3 p 3 = xi {x\ + x 2 3 - d 2 ) + 
+ x 2 (x 3 + x\ — di) + x 3 (x\ + x\ — di), 

p 5 = xi dipi + x 2 d 2 p2 + x 3 d 3 p 3 = Xi di {x\ + x\ - d 2 ) 
+ x 2 d 2 (x% + x\ - di) + x 3 d 3 (xf + x\ - d§), 



2 i 2 | 2 2/2,2 i2 

p e = x 1 p 1 +x 2 p 2 +x 3 p 3 = x 1 (x 2 + x 3 - di 
+ x\ (x 2 +x{- di) + xj (x 2 +x 2 - d 2 ), 

p- = d 2 pi+di p 2 + di p 3 = d\ (x\ + x\ - d\ 
+ d 2 (x 2 +xj- di) + d 2 (x 2 +x 2 - d 2 ), 

= x 2 1 d' 2 p 1 + x\ di p 2 + x\ d'i p 3 = x\ d\ (x\ + x\ 
+ x\ di (x 2 +xj-di)+ x 2 3 d 2 (xj +x 2 -d 2 ). 



(1.12) 
(1.13) 
(1.14) 
(1.15) 
(1.16) 
(1.17) 
(1.18) 



Theorem 1.1. The ideal /pc.sym from (1.10) is finitely generated within the ring 
SymQ[M, L]. Eight polynomials (1.11), (1.12), (1.13), (1.14), (1.15), (1.16), (1.17), 
and (1.18) belong to the ideal ipc_sym o,nd constitute a basis of this ideal. 

The theorem 1.1 was proved in [46]. The ideal ic.sym in (1-8) is similar to the 
ideal Ipc.sym- There is the following theorem describing this ideal. 

Theorem 1.2. The ideal Ic.sym from (1.10) is finitely generated within the ring 
SymQ[M, L\. Seven polynomials (1.12), (1.13), (1.14), (1.15), (1.16), (1.17), and 
(1.18) belong to the ideal /c_sym an d constitute a basis of this ideal. 

The theorem 1.2 can be proved in a way similar to the proof of the theorem 1.1 
in [46]. I do not give the proof of the theorem 1.2 here for the sake of brevity. 

Relying on the theorem 1.2 and using the polynomials (1.12), (1.13), (1-14), 
(1.15), (1.16), (1.17), (1.18), we write the system of seven factor equations 

P2 = 0, p 3 = 0, P4=0, 

P5 = 0, p 6 = 0, p 7 = 0, p 8 = 0. 
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The factor equations (1.19) correspond to the case of Euler cuboids. Similarly, in 
the case of perfect cuboids, relying on the theorem 1.1 and using the polynomials 
given by the formulas (1.11), (1.12), (1.13), (1.14), (1.15), (1.16), (1.17), (1.18), we 
write the following system of eight factor equations: 

Pl=0, P2=0, p 3 = 0, p 4 = 0, 

P5 = 0, Pa = 0, p 7 = 0, p 8 = 0. 

The structure of the polynomials pi, p2, p3, pi, P5, P6, Pi, P8 in (1-H)j (1-12), 
(1.13), (1.14), (1.15), (1.16), (1-17), (1-18) is so that each solution of the equations 

(1.2) is a solution for the equations (1.19). Similarly, each solution of the equations 

(1.3) is a solution for the equations (1.20). The main goal of this paper is to prove 
converse propositions. They are given by the following two theorems. 

Theorem 1.3. Each integer or rational solution of the factor equations (1.19) such 
that xi > ; X2 > 0, X3 > 0, di > 0, e?2 > 0, and d^ > is an integer or rational 
solution for the equations (1-2). 

Theorem 1.4. Each integer or rational solution of the factor equations (1.20) such 
that xi > 0, X2 > 0, X3 > 0, di > 0, e?2 > 0, and d% > is an integer or rational 
solution for the equations (1.3). 



2. The analysis of the factor equations. 

Let's consider the factor equations (1-19) associated with Euler cuboids. Due to 
(1.12), (1.13), (1.14), (1.15), (1.16), (1.17), and (1.18) the factor equations (1.19) 
can be united into a single matrix equation 



1 


1 


di 


d 2 


Xl 


X2 


Xi di 


X 2 d2 






xj 


./'•) 


dl 


dl 


x\ dl 


x\ c?| 













Pi 







Pi 







P3 


















1 

d 3 

X3 

x 3 d 3 

4 

dl 

x\ d| 



In order to study the equations (2.1) we denote through N the transposed matrix 



N = 



1 di xi 
1 d 2 x 2 

1 G?3 X3 



Xi di 
X2 d 2 
x 3 d 3 



r 2 

•' 1 



dl xldl 

X 



2 "2 

2 12 

3 a 3 



(2.2) 



If we have a solution of the equation (2.1) which is not a solution for the initial 
system of cuboid equations (1.2), then the equations (1.2) should not be fulfilled 
simultaneously. Therefore we have the vectorial inequality 



Pi 

P2 
P3 



(2.3) 



ON THE EQUIVALENCE OF CUBOID EQUATIONS 



5 



Applying (2.3) to (2.1), we derive that the columns of the matrix in (2.1) are linearly 
dependent. Then the rows of TV in (2.2) are also linearly dependent, i.e. 



rank AT s$ 2. 



(2.4) 



The condition (2.4) leads to several special cases which are considered below one 
by one. In addition to TV we define the following two matrices: 





1 


di 




1 


Xi 


Ni = 


1 


d 2 


N 2 = 


1 


X 2 




1 


d 3 




1 


X3 



(2.5) 



The matrices N\ and N 2 in (2.5) are submatrices of the matrix TV. 



3. The case rank TV = 1. 

The first column of the matrix (2.2) is nonzero. Therefore rank TV > 0. Now we 
consider the case where rank TV = 1. In this case each column of the matrix TV is 
proportional to its first column. In particular, this yields 



X\ 




1 




di 




1 


X2 


= a ■ 


1 


J 


d 2 


= 13- 


1 






1 




d 3 




1 



The equations (3.1) lead to the equalities 



xi=x 2 =X3 7 di = d 2 = d 3 . (3.2) 

Applying (3.2) to the formulas (1.1), we derive 

P\=p 2 =PZ- (3.3) 

Then we substitute (3.3) into (1.12). As a result we get 

Pi = 3pi = 3p 2 = 3p 3 . (3.4) 

The relationships (3.4) mean that if the equations (1.19) are fulfilled, then in the 
case of rank TV = 1 the equations (1.2) are also fulfilled. 

Theorem 3.1. Each solution of the equations (1.19) corresponding to the case 
rankTV = 1 is a solution for the equations (1.2). 

Theorem 3.2. Each solution of the equations (1.20) corresponding to the case 
rankTV = 1 is a solution for the equations (1.3). 

Due to (1.11) the equation po = in (1.3) coincides with the equation pi = in 
(1.20). For this reason the theorem 3.2 is immediate from the theorem 3.1. 

4. The case rankTVi = 2 and rankTV 2 = 1. 

The condition rankTV 2 = 1 for the matrix N 2 in (2.5) means that the third 
column of the matrix (2.2) is proportional to the first column of this matrix. The 
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condition rank Ni = 2 for the matrix Ni in (2.5) means that the first and the second 
columns of the matrix (2.2) are linearly independent. Other columns are expressed 
as linear combinations of these two columns. As a result we can write 



(4.1) 



It is easy to see that the conditions (4.1) are sufficient for to provide the condition 
rank./V = 2, which is in agreement with (2.4). 

The second equality in (4.1) is very important. It means that d\, di-, and d 3 , are 
roots of the following quadratic equation: 



Xi 




1 


X2 


= a ■ 


1 


X 3 




1 



dl 




di 




1 


dl 


= P- 


d 2 


+ 7' 


1 


dl 




d 3 




1 



(d-7 = 0. 



(4.2) 



The quadratic equation (4.2) has at most two roots. Let's denote them si and S2- 
Then we have the following subcases derived from rank N\ = 2 and rank N% = 1: 



(4.3) 



The numbers Si and S2 in the formulas (4.3) are arbitrary two numbers not coin- 
ciding with each other: si ^ S2- They are integer numbers in the case of integer 
solutions and they are rational numbers in the case of rational solutions. 

The three cases in (4.3) are similar to each other. Without loss of generality we 
can consider only one of them, e. g. the first one. Then from (4.3) we derive 



di 




Sl 




di 




Sl 


d 2 




Sl 


1 


d 2 




S2 


d 3 




s 2 




d 3 




Sl 



di 




S2 


d 2 




Sl 


d 3 




Sl 



(si - s 2 ) 



(S2 ~ Sl) 



1 




Sl 




1 


1 




Sl 


- s 2 ■ 


1 







32 




1 







Sl 




1 







Sl 


- Sl ■ 


1 


1 




S2 




1 



(4.4) 



(4.5) 



Due to the relationships (4.4) and (4.5) the matrix equation (2.1) reduces to 



1 


1 







Pi 













1 




Pi 
p 3 








(4.6) 



The matrix equality (4.6) means that instead of the seven equations (1.19) we 
have two equations pi + P2 =0 and p 3 = 0. Substituting xi = X2 = £3 = a, 
di = c?2 = si, and d 3 = S2 into these two equations, we derive 



s\ - 2 a 2 = 0, 



s 2 -2 a 2 



0. 



(4.7) 



The equations (4.7) can be written in the following way: 



\si\ = V2\a\, 



\s 2 \ = V2\a\. 



(4.8) 
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Now it is easy to see that the equations (4.8) can be satisfied by three integer or 
rational numbers si, s 2 , and a if and only if all of them are zero. Substituting 
si = S2 = a = into (4.1) and (4.3), we get 



xi 



X2 = X 3 = 0, 



di 



0. 



(4.9) 



The equalities (4.9) contradict the condition rankiVi = 2 for the matrix N\ in (2.5). 
This contradiction yields the following two theorems. 

Theorem 4.1. The factor equations (1.19), as well as the original equations (1.2), 
have no integer or rational solution in the case o/rankiVi = 2 and rankiV2 = 1. 

Theorem 4.2. The factor equations (1.20), as well as the original equations (1.3), 
have no integer or rational solution in the case o/rankiVi = 2 and rankiVa = 1. 

5. The case rankiVi = 1 and rank7V 2 = 2. 

The condition rankiVi = 1 for the matrix TVi in (2.5) means that the second 
column of the matrix (2.2) is proportional to the first column of this matrix. The 
condition rankiV 2 = 2 for the matrix N2 in (2.5) means that the first and the 
third columns of the matrix (2.2) are linearly independent. Other columns are 
expressed as linear combinations of these two columns. As a result we can write 
the relationships similar to the relationships (4.1): 



di 




1 


d 2 


= 5- 


1 


d 3 




1 



i 




Xi 




1 




= e ■ 


x-i 


+ c- 


1 


T 2 




£3 




1 



(5.1) 



The conditions (5.1) are sufficient for to provide the condition rankiV = 2. 

Like in the case of (4.1), the second condition (5.1) mean that xi, X2, and £3 are 
roots of the quadratic equation similar to (4.2): 



C = o. 



(5.2) 



The quadratic equation (5.2) has at most two roots. Let's denote them r\ and ti. 
Then we have the following subcases derived from rankiVi = 1 and rankA^ = 2: 



Xl 




n 




Xl 




n 


X2 




n 


1 


X2 






%3 




T2 




X3 




n 



Xl 




T2 


X2 




n 


%3 




n 



(5.3) 



The numbers r\ and r2 in the formulas (5.3) are arbitrary two integer or rational 
numbers not coinciding with each other: r\ 7^ 

The three cases in (5.3) are similar to each other. Without loss of generality we 
can consider only one of them, e.g. the first one. Then from (5.3) we derive 





1 




n 




1 


(n - r 2 ) • 


1 




n 


- r 2 ■ 


1 









T2 




1 



(5.4) 
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n 




l 


(r-2 - n) ■ 







n 


- n ■ 


l 




1 




T2 




l 



(5.5) 



Due to the relationships (5.4) and (5.5) the matrix equation (2.1) reduces to 



1 


1 







Pi 













1 




P2 















P3 







(5.6) 



The matrix equality (5.6) means that instead of the seven equations (1-19) we have 
two equations p i +p 2 = andpz = 0. Substituting d\ = d% = ds = 5, x% = Xi = r%, 
and xz = f2 into these two equations, we derive 



rj+rj-5 2 =0, 2rj-S 2 = 0. 

The second equation (5.7) can be written in the following form: 

|*| = V2|n|. 



(5.7) 



(5. 



The equation (5.8) can be satisfied by two integer or rational numbers n and S if 
and only if both of them are zero. Substituting n = 6 = into the first equation 
(5.7), we get = 0. Substituting n = r 2 = S = into (5.1) and (5.3), we get 



Xl 



X2 = X3 = 0, 



0. 



(5.9) 



The equalities (5.9) contradict the condition rank = 2 for the matrix N2 in (2.5). 
This contradiction yields the following two theorems. 

Theorem 5.1. The factor equations (1.19), as well as the original equations (1-2), 
have no integer or rational solution in the case of rank Ni = 1 and rankiV2 = 2. 

Theorem 5.2. The factor equations (1.20), as well as the original equations (1-3), 
have no integer or rational solution in the case of rank iVi = 1 and rank7V2 = 2. 

6. The case rankTVi = 2 and rank7V 2 = 2. 

In this case the columns of both matrices Ni and A^ 2 in (2.5) are linearly in- 
dependent. Hence each column of the matrix (2.2) can be expressed as a linear 
combination of the first and the second columns of this matrix or as a linear com- 
bination of the first and the third columns of this matrix. In particular, we have 



(6.1) 



The relationships (6.1) mean that x\ , X2 , ^3 and d\ , c?2 , ^3 are roots of two quadratic 
equations coinciding with (5.2) and (4.2) respectively. As a result we distinguish 
three subcases (4.3) with s\ ^ s 2 and three subcases (5.3) with n ^ ?'2. The first 
subcase (4.3) should be paired with the first subcase (5.3), the second subcase (4.3) 
should be paired with the second subcase (5.3), and the third subcase (4.3) should 



dl 




d! 




1 




A 




Xl 




1 


4 


= P- 


d 2 


+ 7' 


1 




x 2 


= e ■ 


X2 


+ c- 


1 


d 2 




d 3 




1 




T 2 

x 3 




X3 




1 
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be paired with the third subcase (5.3). Otherwise we would have rank TV ^ 3, which 
contradicts the condition (2.4). 

Due to the pairing of subcases we have three subcases instead of nine ones, which 
are a priori possible. These three subcases are similar to each other. Therefore 
without loss of generality we can consider only one subcase, e. g. the following one: 



(6.2) 



Here si ^ S2 and r\ ^ ri. The relationships (6.2) lead to the relationships (4.4), 

(4.5) , (5.4), (5.5) and then to the equations (4.6) and (5.6). The matrix equations 

(4.6) and (5.6) mean that instead of the seven equations (1.19) we have two equa- 
tions pi + P2 = and P3 = 0. Substituting d\ = c?2 = si, d% = S2, x± = x<i = n, 
and X3 = T2 into these two equations, we derive 



d x 




Sl 




Xl 




ri 


d 2 




81 


1 


X2 




ri 


d 3 




S2 




X3 




T2 



r\+rl- s\ = 0, 



2r? -s 2 = 0. 



(6.3) 



The second equation (6.3) can be written in the following way: 

V2\n\ = \s 2 \. (6.4) 

The equation (6.4) can be satisfied by two integer or rational numbers n and s 2 
if and only if both of them are zero. Substituting n = S2 = into (6.3), we get 
1 7^2 1 = | s 1 1 = 6. Then the equalities (6.2) are written as 



(6.5) 



The equalities (6.5) lead to the equalities xi = X2 — and d 3 = 0. The latter ones 
contradict the inequalities in the theorems 1.3 and 1.4. Therefore we can conclude 
this section with the following two theorems. 

Theorem 6.1. The factor equations (1.19), as well as the original equations (1-2), 
have no integer or rational solutions such that X\ > 0, X2 > 0, X3 > 0, d\ > 0, 
c?2 > 0, and d 3 > in the case of rank N± = 2 and rankA^ = 2. 

Theorem 6.2. The factor equations (1.20), as well as the original equations (1.3), 
have no integer or rational solutions such that x\ > 0, X2 > 0, X3 > 0, di > 0. 
c?2 > 0, and c?3 > in the case o/rankiVi = 2 and rankA^ = 2. 



di 




e 




Xl 







d 2 


= ± 


e 


1 


X2 


= ± 





d 3 









X?, 
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7. The ultimate result and conclusions. 

The four cases considered in sections 3, 4, 5, and 6 exhaust all options compatible 
with the inequality (2.4). For this reason the theorems 1.3 and 1.4 follow from the 
theorems 3.1, 4.1, 5.1, 6.1 and the theorems 3.2, 4.2, 5.2, 6.2 respectively. The 
theorems 1.3 and 1.4 constitute the main result of this paper. The theorem 1.4 
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means that the factor equations (1.20) arc equally admissible for seeking perfect 
cuboids or for proving their non-existence as the original equations (1.3). As for 
the factor equations (1.19), due to the theorem 1.3 they are equally admissible for 
selecting Euler cuboids as the original equations (1.2). 
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